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Abstract 

We present results for the 3 P2- 3 -^2 pairing gap in neutron matter with sev- 
eral realistic nucleon-nucleon potentials, in particular with recent, phase-shift 
equivalent potentials. We find that their predictions for the gap cannot be 
trusted at densities above p « 1.7po, where po is the saturation density for 
symmetric nuclear matter. In order to make predictions above that density, 
potential models which fit the nucleon-nucleon phase shifts up to about 1 GeV 
are required. 
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I. INTRODUCTION 



The presence of neutron superfluidity in the crust and the inner part of neutron stars 
is one of the features that are considered well established in the physics of these compact 
stellar objects. At low density, and therefore in the outer part of a neutron star, the neutron 
superfluidity should be mainly in the 1 S'o channel. At higher density, the nuclei in the crust 
dissolve, and one expects a region consisting of a quantum liquid of neutrons and protons 
in beta equilibrium. The proton contaminant should be superfluid in the 1 Sq channel, 
while neutron superfluidity is expected to occur mainly in the coupled 3 P2~ 3 F 2 two-neutron 
channel. In the core of the star any superfluid phase should finally disappear. 

The presence of two different superfluid regimes is suggested by the known trend of the 
nucleon-nucleon (NN) phase shifts in each scattering channel. In both the 1 5'o and 3 P 2 - 3 F 2 
channels the phase shifts indicate that the NN interaction is attractive. In particular for the 
1 S'o channel, the occurrence of the well known virtual state in the neutron-neutron channel 
strongly suggests the possibility of a pairing condensate at low density, while for the 3 P 2 - 3 F 2 
channel the interaction becomes strongly attractive only at higher energy, which therefore 
suggests a possible pairing condensate in this channel at higher densities. In recent years 
the BCS gap equation has actually been solved with realistic interactions, and the results 
confirm these expectations. 

The 1 Sq neutron superfluid is relevant for phenomena that can occur in the inner crust of 
neutron stars, like the formation of glitches, which seem to be related to vortex pinning of the 
superfluid phase in the solid crust QXJ . The results of different groups are in close agreement 
on the 1 5'o pairing gap values and on its density dependence, which shows a peak value of 
about 3 MeV at a Fermi momentum close to kp m 0.8 fm _1 0-0. All these calculations 
adopt the bare NN interaction as the pairing force, and it has been pointed out that the 
screening by the medium of the interaction could strongly reduce the pairing strength in 
this channel 0-0] . However, the issue of the many-body calculation of the pairing effective 
interaction is a complex one and still far from a satisfactory solution. 

The precise knowledge of the 3 P2- 3 F2 pairing gap is of paramount relevance for, e.g., the 
cooling of neutron stars, and different values correspond to drastically different scenarios 
for the cooling process ||. Unfortunately, only few and partly contradictory calculations 
of this quantity exist in the literature, even at the level of the bare NN interaction PHTB" 



However, when comparing the results, one should note that the NN potentials used in these 
calculations are not phase-shift equivalent, i.e., they do not predict exactly the same NN 
phase shifts. Furthermore, for the interactions used in Refs. pj-[12[| the predicted phase shifts 
do not agree accurately with modern phase shift analyses, and the fit of the NN data has 
typically xVdatum w 3. During the last years, progress has been made not only in the 
accuracy and the consistency of the phase-shift analysis, but also in the fit of realistic NN 
potentials to these data. As a result, several new NN potentials have been constructed which 
fit the world data for pp and np scattering below 350 MeV with high precision. Potentials 
like the recent Argonne V 18 the CD-Bonn or the new Nijmegen potentials yield 



a % 2 /datum of about 1 and may be called phase-shift equivalent. 

Our aim in this paper is to compare the predictions of the new potentials for the 3 P 2 - 3 F 2 
gap in neutron matter. We will also, for the sake of completeness, include results with some 
of the "old" interactions, namely the Paris [17] , Argonne V14 |1| , and Bonn B |19| potentials. 
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The main focus will, however, be on the new, phase-shift equivalent potentials, and whether 
the improved accuracy in the fits to the NN scattering data leads to better agreement in the 
predictions for the 3 P 2 - 3 F 2 energy gap. If differences are found, we try to trace them back 
to features of the NN potentials. To be able to do so, we will keep the many-body formalism 
as simple as possible. First of all, we will use the bare NN interaction as kernel in the 
gap equations, and thus neglect higher-order contributions from, e.g., medium polarization 
effects. The in-medium single-particle energies will be calculated in the Brueckner-Hartree- 
Fock (BHF) approximation, but we will also use free single-particle energies, because this 
makes the comparison of the results with the various potentials more transparent, since any 
differences are then solely due to differences in the 3 P 2 - 3 F 2 wave of the potentials. We think 
it is useful to try to understand the results at the simplest level of many-body theory before 
proceeding to include more complicated effects in the description of 3 P 2 - 3 F 2 pairing. As we 
will demonstrate, progress in the construction of NN interactions is necessary before the 
3 P 2 - 3 F 2 energy gap can be calculated reliably from microscopic many-body theory. 

This work falls in six sections. The equations for solving the pairing gap are briefly 
reviewed in the next section, while in Section |T| we discuss the reliability of various numerical 
approaches to the solution of the pairing gap. Features of the various nucleon-nucleon 
interaction models employed are presented in Section [TV], while our results for the pairing 
gap with these potentials are discussed in Section [V|. Finally, we summarize our findings in 



Section VI 



II. GAP EQUATION FOR THE 6 P 2 -F 2 CHANNEL 

The gap equation for pairing in non-isotropic partial waves is in general more complex 
than in the simplest s-wave case, in particular in neutron and nuclear matter, where the 
tensor interaction can couple two different partial waves [PH , ^0[| . This is indeed the situation 
for the 3 P 2 - 3 F 2 neutron channel. In order to achieve a simplified, yet accurate, numerical 
treatment, we use in this work the angle average approximation explained in this section. 

For the sake of a clear presentation, we disregard for the moment the spin degrees of 



freedom and the tensor interaction. Then, we start with the Gorkov equations [pi|| , which 
involve the propagator G(k,u), the anomalous propagator F(k,u), and the gap function 
A(fc): 



oo-e(k) -A(fc)\ / G 
-At(Jfe) iu + e(k) ) (Ft 



> (1) 



where e(k) = e(k) — /i, /i being the chemical potential and e(k) the single-particle spectrum. 
The quasi-particle energy E(k) is the solution of the corresponding secular equation and is 
given by 

£(fc) 2 = e (fc) 2 +|A(fc)| 2 . (2) 
The anisotropic gap function A(k) is to be determined from the gap equation 

A(*H-X>IW^). (3) 
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The angle-dependent energy denominator in this equation prevents a straightforward sepa- 
ration into the different partial wave components by expanding the potential, 

(k\V\k') = 4vr^(2/+ l)P,(fc • fc'MMO , (4) 

i 

and the gap function, 



4-7T 

^(k)=Y /] J ¥TI Y lm (k)A lm (k). (5) 

However, after performing an angle average approximation for the gap in the quasi-particle 
energy, 

|A(fc)| 2 -> D(k) 2 = ^~jdk |A(fc)| 2 = £ ^Tll^WI 2 > ( 6 ) 

the kernels of the coupled integral equations become isotropic, and one can see that the dif- 
ferent m-components become uncoupled and all equal. One obtains the following equations 
for the partial wave components of the gap function: 

A, w= -ir , mm a, w . ( 7) 

Note that there is no dependence on the quantum number m in these equations, however, 
they still couple the components of the gap function with different I ( 1 S , , 3 Po, 3 Pi, 3 P 2 , 
1 D 2 , 3 F 2 , etc. in neutron matter) via the energy denominator. Fortunately, in practice the 
different components Vi of the potential act mainly in non-overlapping intervals in density, 
and therefore also this coupling can usually be disregarded. 

The addition of spin degrees of freedom and of the tensor force does not change the 



picture qualitatively, and is explained in detail in Refs. |Tl|^0j. The only modification is the 
introduction of an additional 2x2 matrix structure due to the tensor coupling of the 3 P 2 
and 3 F 2 channels: 

E(k) 2 = [e(k) — e(k F )} 2 + D(k) 2 , (8b) 
D(k) 2 = A^k) 2 + A 3 (k) 2 . (8c) 

Here e(k) = k 2 /2m + U(k) are the single-particle energies, as obtained from a Brueckner- 
Hartree-Fock calculation, where U(k) is the single-particle potential, calculated within the 
"continuous choice" scheme |2"2fl. The quantities 



V w (k,k')= drr 2 ]v (k'r)V ni (r) 3l (kr) (9) 



o 



are the matrix elements of the bare interaction in the different coupled channels (T — 1; S = 
1; J = 2; 1,1' = 1,3). 

It has been shown that the angle average approximation is an excellent approximation 
to the true solution that involves a gap function with ten components |Tl|JT3]] , as long as one 



is only interested in the average value of the gap at the Fermi surface, = D(k F ), and 
not the angular dependence of the gap functions Ai(fe) and A 3 (fc). 
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III. NUMERICAL SOLUTION 



The solution of the system of equations (|8|) is numerically not trivial, especially if the 
gap turns out to be much smaller than the Fermi energy. This is because of the well known 
logarithmic singularity of the BCS equation in the limit of zero pairing gap. In order to 
control more closely the numerical accuracy, we used in fact three different methods: 

One method is similar to the one described in Ref. 0. We first obtain a separable form 
of the interaction. Since we need a high accuracy, we directly diagonalize the interaction 
Vui(k, k'), taken in a discrete grid of momenta {fcj}, and then we choose the first n eigenvalues 
X m with the largest moduli, and the corresponding eigenvectors v m . One can then write 

n 

V lv {k h k-) « v m{k i )\ m v m {k j ) . (10) 

m=l 

The gap function can then also be expanded in the same eigenvectors, and the original equa- 
tions reduce to a set of 2n algebraic equations. The latter can be solved for the coefficients 
of the expansion by iteration, following the scheme described in Ref. || . The rank n of the 
separable form is increased until a high degree of convergence is reached. One advantage 
of the method is the possibility of using a very fine momentum grid, since the algebraic 
equations are obtained by numerical integrations, for which extremely accurate interpola- 
tion methods can be used. In general, the grid points must be particularly dense in the 
interval around the Fermi momentum, since there the kernel displays an extremely narrow 
peak due to the small value of the pairing gap. Furthermore, in general, the convergence 
in the rank n is fast enough, and therefore the number of coupled equations is never very 
large. However, the accuracy in the diagonalization procedure is decreasing with the rank 
of the matrix and it is difficult to have a precise estimate of the error. 

In the second method [E3J one starts by solving the gap equation for the case of a constant 



pairing gap A in the denominator. In a discrete momentum grid, this is equivalent to an 
eigenvalue problem, namely to find the value of A for which the kernel of the gap equation 
has eigenvalue one. The corresponding eigenvector is a first estimate of the gap function, 
with the normalization A(k F ) = A. It is then inserted in the kernel to solve for the next 
estimate of A. In practice this method converges extremely fast (after a few iterations) to 
the final solution. The advantage of the method is that the original interaction is used, 
without resorting to a separable form. 

The third method is to solve the coupled 3 P 2 - 3 F 2 gap equations straightforwardly by 
iteration, starting from some suitable initial approximation to the functions Ai(fc) and 
A 3 (/c). Also in this method, the interaction is used in its original form. If the interaction 
has a strong repulsive core, as is the case in the 1 5'o channel, this method can be difficult 
or even impossible to implement. However, the 3 P 2 - 3 F 2 interaction is relatively weak, and 
the iteration scheme works well in this channel, provided that a fine momentum grid is used 
around the Fermi momentum. Details of the numerical implementation of this method are 



given in Ref. [|12 



The comparison of the results obtained with the three methods was quite rewarding. The 
numerical values of the gap functions were in excellent agreement and hardly distinguishable 
in all figures presented here. Therefore, in discussing the results we will not specify the 
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method by which they were obtained. We believe that the agreement between the three 
methods gives enough confidence in the numerical precision of the results. 



IV. THE NN INTERACTIONS 

Before discussing the solutions of the coupled 3 P2- 3 F 2 gap equations, we give a short 
description of the models for the NN interaction employed in this paper. 



The older models, Paris, Argonne V14, and Bonn B are described in detail in Refs. fl7|-|19 
They all have a x 2 /datum in the range 2-3. The Argonne Vu potential is a non-relativistic, 
purely local potential. The Paris potential incorporates explicit 7T-, 2tt-, and u;-exchange. 
For the short-range part a phenomenological approach is used. The final potential is param- 
eterized in terms of local Yukawa functions. The Bonn B potential is a one-boson-exchange 
(OBE) interaction, defined by the parameters of Table A.l of Ref. |L9|. 



The "phase-shift equivalent" potentials we will employ here are the recent models of 
the Nijmegen group |16|, the Argonne Vis |H[ potential and the charge- dependent Bonn 



potential (CD-Bonn) [JTJJ. In 1993, the Nijmegen group presented a phase-shift analysis of 
all proton-proton and neutron-proton scattering data below 350 MeV with a xVdatum of 
0.99 for 4301 data entries j23j. Fitted to this phase-shift analysis, the CD-Bonn potential 
has a xVdatum of 1.03 and the same is true for the Nijm-I and Nijm-II potentials of the 
Nijmegen group [16|]. The Argonne Vis potential has a x /datum of 1.09. 

All these models are charge-dependent. Argonne Vis and Nijm-II are non-relativistic 
potential models defined in terms of local functions, which are attached to various (non- 
relativistic) operators constructed from the spin, isospin and angular momentum operators of 
the interacting pair of nucleons. Such approaches to the NN potential have traditionally been 
quite popular since they are numerically easy to use in configuration space calculations. The 
Nijm-I model is similar to the Nijm-II model, but it includes also a momentum dependent 
term, see Eq. (13) of Ref. ||16|| , which may be interpreted as a non-local contribution to the 



central force. The CD-Bonn potential is based on the relativistic meson-exchange model of 



Ref. [|19j which is non-local and cannot be described correctly in terms of local potential 
functions. Instead, it is represented most conveniently in terms of partial waves. 

Thus, the mathematical structure of the modern potentials is quite different, although 
they all predict almost identical phase shifts within their range of validity. This means that 
even though the potentials by construction give the same results on-shell, their behavior 
off the energy shell may be quite different. The implications of these differences for the 
symmetry energy of nuclear matter were discussed in Ref. [^5] . 



In order to illustrate the statements made above, and for a better understanding of the 
forthcoming results for the pairing gaps, we show in Fig. |l] the predictions of the various 
potentials for the phase shifts in the 3 P2 (T = 1) channel. They have been calculated 
by solving the Lippmann-Schwinger equation as explained in Ref. [p^j . The figure shows 
predictions up to = 1.1 GeV, but clearly scattering energies above E^b = 350 MeV 
amount to uncontrolled extrapolations beyond the intended range of validity of the potential 
models, that have been fitted to scattering data below 350 MeV only. The plot displays also 
a scale of equivalent Fermi momenta according to the relation E^ = (2/ci?) 2 /2m in order to 
facilitate the comparison with the pairing gaps presented later. The reader can see that a lab 
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energy of 350 MeV corresponds roughly to a Fermi momentum kp = 2.0 fm~ . Therefore, 
calculations of the 3 P 2 - 3 F 2 energy gap at densities above kp = 2.0 fm _1 will inevitably 
involve extrapolating the potential models. 

In the same figure we also show the empirical pp phase shifts obtained by Arndt et al. 
in a recent phase shift analysis [p7 |. Some differences between this phase shifts analysis 
and the phase shifts calculated with the potentials could be present in the figure, even 
below 350 MeV, because the potentials are not fitted to the analysis of Arndt et al.. The 



modern potentials fits the Nijmegen database |2J], the older ones fit different analyses made 



in the 70's and 80's. Nevertheless, the four modern potentials considered here fit also Arndt 
analysis below 350 MeV with high accuracy, while the old potentials (in particular the V14) 
overshoot the empirical values already at lower scattering energies, due to the fact that they 
have a higher x 2 /datum than the new models. 

In any case, above E^ = 350 MeV (corresponding to kp « 2.0 fin -1 ) sizeable differences 
show up in the predictions of all potentials. The Nijm-II potential fits the phase shifts up 
to about 600 MeV rather well, but after that it severely overestimates them. This in turn 
means that the high-momentum components of the 3 P 2 interaction will be too attractive. 
Nijm-I does fairly well up to about 500 MeV, from 500 to 700 MeV it underpredicts the 
phase shifts, while at energies above 700 MeV the results are too high. The CD-Bonn 
potential gives a similar behavior, but falls faster towards zero at high energies than Nijm-I 
and II. Argonne Vi$ gives 3 P 2 phase shifts below the empirical ones over the whole range 
E\ a h — 400-1000 MeV. The old potentials display similar variations, being generally too 
repulsive with Paris the most repulsive of all potentials, followed by Bonn B and Argonne 
V14. In this paper we will further on focus on the new, phase-shift equivalent potentials, 
since they are better fitted to modern scattering data. In summary, all potentials give phase 
shifts which are too attractive above ph 700-1000 MeV, and all except Nijm-II are too 
repulsive between « 350 MeV and « 700-1000 MeV. 



V. RESULTS 

Before presenting results for the energy gap, we point out some features of the gap equa- 
tions which make the trend of the results understandable. In order to make the connection 
to the NN interaction as transparent as possible, we start by discussing the case where the 
single-particle energies are given by their values in free space, e(k) = k 2 /2m. 

In Fig. |2| we show, for the Nijm-I potential and various values of kp, the function 
k 2 Ai(k)/E(k) involved in the 3 P 2 component of the gap equations, normalized to unity 
at k = kp. The behavior of this function was found to be the same for all potentials. No- 
tice that this function is very strongly peaked around k = kp, implying that the diagonal 
matrix element of the potential at k — kp gives the most important contribution to Ai(kp) 
and As(kp). Also, this figure makes it clear why some care in choosing momentum mesh 
points for the numerical integrations is needed. The function k 2 A 3 (k)/E(k) shows a similar, 
strongly peaked behavior, and thus the gap is largely determined by the matrix elements 
V xl {k F , kp), V 13 (k F ,k F ) and V 33 (kp,k F ). 

To exemplify this, we have therefore plotted in Fig. |^ the matrix elements for V\\{kp, kp) 
and V 33 (kp, kp) as functions of kp for the various modern potentials used in this work. Up 
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to kp « 2.0 fm _1 the matrix elements are very similar, but after this point they deviate 
from each other, in line with the phase shift predictions shown in Fig. [TJ: In the 3 P 2 and 3 P 2 
waves, the Vi§ potential is the most repulsive, followed by the CD-Bonn and the Nijm-I and 
Nijm-II potentials in that order. Similar conclusions can be reached for the coupled 3 P 2 - 3 F 2 
channel. 



A. Pairing gaps 

Fig. |] contains a comprehensive collection of our results for the pairing gaps with the 
different potentials. We start with the top part of the figure, which displays the results 
calculated with free single-particle energies. Differences between the results are therefore 
solely due to differences in the 3 P 2 - 3 F 2 matrix elements of the potentials. The plot shows 
results obtained with the old as well as with the modern potentials. The results (with 
the notable exception of the Argonne V14Q, which predicts also substantially different 3 P 2 
phase shifts (see Fig. []]), are in good agreement at densities below kp 2.0 fm _1 , but 
differ significantly at higher densities. This is in accordance with the fact that the diagonal 
matrix elements of the potentials are very similar below kp rs 2.0 fm -1 , corresponding to a 
laboratory energy for free NN scattering of E^ rs 350 MeV. This indicates that within this 
range the good fit of the potentials to scattering data below 350 MeV makes the ambiguities 
in the results for the energy gap quite small, since, to a first approximation, see the discussion 
below, the pairing gap can be derived in terms of the phase shifts only. 

However, we wish to calculate the gap also at densities above kp = 2.0 fm _1 . Then we 
need the various potentials at higher energies, outside of the range where they are fitted 
to scattering data. Thus there is no guarantee that the results will be independent of the 
model chosen, and in fact the figure shows that there are considerable differences between 
their predictions at high densities, following precisely the trend observed in the phase shift 
predictions: The Argonne Vi§ is the most repulsive of the modern potentials, followed by 
the CD-Bonn and Nijmegen I and II. Most remarkable are the results obtained with Nijm- 
II: we find that the predicted gap continues to rise unrealistically even at kp ~ 3.5 fm _1 , 
where the purely nucleonic description of matter surely breaks down. From Table |, which 
contains a compilation of gaps for the various potentias, one sees that the improved fit of 
the new potentials to scattering data leads to better agreement in their predictions for the 
gap. Thus, the fact that these potentials have been fit with high precision to the same set of 
scattering data eliminates some of the ambiguities, and allows one to compare interactions 
in a way not possible with earlier models. 

Since the potentials fail to reproduce the measured phase shifts beyond = 350 MeV, 
the predictions for the 3 P 2 - 3 F 2 energy gap in neutron matter cannot be trusted above kp ~ 
2.0 fin -1 . Therefore, the behavior of the 3 P 2 - 3 P 2 energy gap at high densities should be 
considered as unknown, and cannot be obtained until potential models which fit the phase 



1 In a previous paper [ 10 1 one of the authors (M.B.) has claimed much higher values for the gap 
with the Argonne V14. It has been checked that this was due both to a non accurate separable 
representation of the NN potential and to a bug in the computer program for this channel. 
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shifts in the inelastic region above E\ a ^ = 350 MeV are constructed. These potential models 
need the flexibility to include both the flat structure in the phase shifts above 600 MeV, due 
to the NN — > NA channel, as well as the rapid decrease to zero at E^b ~ 1100 MeV. 

We proceed now to the middle part of Fig. |], where the results for the energy gap 
using BHF single-particle energies are shown. For details on the BHF calculations, see, e.g., 



Ref. p2| . From this figure, two trends are apparent: First, the reduction of the in-medium 
nucleon mass leads to a sizeable reduction of the 3 P2- 3 P2 energy gap, as observed in earlier 
calculations [PHl2||. Secondly, the new NN interactions give again similar results at low 
densities, while beyond k F « 2.0 fm -1 the gaps differ, as in the case with free single-particle 
energies. 

The single-particle energies at moderate densities obtained from the new potentials are 
rather similar, particularly in the important region near hp- This is illustrated by a plot, 
Fig. [|, of the neutron effective mass, 

/ \ _1 

77T 





dU 






dk 





(11) 



as a function of density. Up to k F rs 2.0 fm -1 all results agree very closely, but beyond that 
point the predictions diverge in the same manner as observed for the phase shift predictions. 
The differences of the BHF gaps at densities slightly above kp ~ 2.0 fm~ are therefore 
mostly due to the differences in the 3 P2- 3 P2 waves of the potentials, but at higher densities 
the differences between the gap are enhanced by differences in the single-particle potentials. 
The reader should bear in mind that the single-particle energies contain contributions from 
partial waves up to I < 10. The largest differences arise however from contributions from 
the l So and 3 P2- 3 F 2 partial waves, see also the discussion in Ref. [25]. An extreme case 
is again the gap obtained with Nijm-II. It is caused by the very attractive 3 P2 matrix 
elements, amplified by the fact that the effective mass starts to increase at densities above 
kp fa 2.5 fm -1 with this potential. 

Finally, in the lower panel of Fig. we illustrate the effect of different approximation 
schemes with an individual NN potential (CD-Bonn), namely we compare the energy gaps 
obtained with the free single-particle spectrum, the BHF spectrum, and an effective mass 
approximation, 

= (12) 

where m* is given in Eq. (|TT|). In addition, also the gap in the uncoupled 3 P2 channel, i.e., 
neglecting the tensor coupling, is shown. 

It becomes clear from Fig. ^ that the BHF spectrum forces a reduction of the gap by about 
a factor 2-3. However, an effective mass aproximation should not be used when calculating 
the gap, because details of the single-particle spectrum around the Fermi momentum are 
important in order to obtain a correct value. The single-particle energies in the effective 
mass approximation are too steep near kp. We also emphasize that it is important to solve 
the coupled 3 P 2 - 3 P2 gap equations. By turning off the 3 P 2 - 3 P2 and 3 P 2 channels, one obtains 
a 3 P2 gap that is considerably lower than the 3 P2- 3 F 2 one. The reduction varies with the 
potential, due to different strengths of the tensor force. For more detailed discussions of the 
importance of the tensor force, the reader is referred to Refs. P JTT|JT^ ]. 



9 



B. Hints from the 3 P2 phase shifts 



The first calculation of the 3 P 2 gap in neutron matter was carried out by Hoffberg et 
al. |28[] in 1970. They used the weak-coupling expression for the energy gap to express it in 
terms of the 3 P 2 phase shifts available at that time, and obtained a maximum gap of around 
1 MeV at kp ~ 2.3 fm . Since all interactions considered in the present paper are fitted in 
the energy range 0-350 MeV, it would be interesting to use the recent phase shift analysis 
by Arndt et al. [^7]] to get some hints on the behavior of the energy gap at higher densities. 
The phase shifts determine the interaction only on the energy shell, so to go from these 
"experimental" points to the energy gaps, we must make some rather strong assumptions. 

First of all, we switch off the interaction in the 3 F 2 and 3 P 2 - 3 F 2 channels and consider 
pure 3 P 2 pairing. We are then left with only one gap equation to solve, and when we use 
the angle average approximation it is identical in form to the equation for 1 5'o pairing: 



A 1 (fc) = ~jf ) dk'k' 2 V n (k,k')^l. (13) 

In a recent paper J3J two of the authors derived an expression for the 1 So gap in neutron and 
nuclear matter in terms of the phase shifts in this partial wave. This was possible because 
the interaction in this channel is to a good approximation rank-one separable at low energies 
due to the 1 Sq two-nucleon virtual state |||29|]. No resonance or virtual state exists in the 3 P 2 



channel, but we will nevertheless approximate the interaction in this channel by a rank-one 
separable form, 

Vn(k,k') = Xv(k)v(k') , (14) 
where A is a constant. The interaction can then be expressed in terms of the phase shifts as 



Xv \k) = _^M e - Q W , (15) 
k 



where a(k) is given by a principle value integral 

7T J-oo k — k 

and the phase shifts 5(k) are extended to negative momenta through 5(—k) = —5(k). This 
prescription works only for a purely attractive or purely repulsive interaction. The 3 P 2 
phase shifts change sign at ~ 1100 MeV, and thus the interaction goes from attractive 
to repulsive at this energy. We therefore cut the integral in Eq. (|i"6|) at k « 3.6 fm -1 , which 
corresponds to £?i a b ~ 1100 MeV. For a rank-one separable interaction, the solution of 
Eq. (|i~3l ) is given by Apv(k), where is the gap at the Fermi momentum found by solving 



Using phase shifts from the analysis of Arndt et al. p7| , |3Tl| , we constructed an interaction 
for the 3 P 2 channel according to the prescription above, and then proceeded to solve Eq. ([H 
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for Ap. The results are shown in Fig. |6|. For comparison we also display the results of the 
following calculation for the various potentials: we took the phase shifts at energies up to 
1100 MeV computed earlier and shown in Fig. [I]. From these we constructed a rank-one 
separable approximation to the 3 P 2 wave of the various potentials, as described above, and 
then used this to solve the gap equation. As such, we have a as close as possible link with 
the calculation based solely on the phase shifts of Arndt et al. |27| , |3T| . This allows us in 
turn to see directly the consequences of the failure of the potentials to fit the high-energy 
3 P 2 phase shifts. When looking at Fig. ^ and reading the following discussion, one should 
bear in mind that the gap has an exponential dependence on the interaction, so quite small 
differences in the matrix elements of the interaction can be translated into large differences 
in the energy gap. But this also makes the gap a good quantity to use when comparing 
interactions, as any difference is magnified. 

Although the approximation made here should not be taken too seriously, the results 
indicate some important conclusions about the 3 P 2 waves of the recent nucleon-nucleon 
interactions. All seem to have about the right amount of attraction at densities below 
kp « 2.0 fm~ . Between hp ~ 2.0 fm~ and kp « 3.0 fm~ all interactions except Nijm-II 
seem to be a bit too repulsive. Above kp « 3.0, Argonne Vis is probably too repulsive, while 
Nijm I and II are most certainly too attractive, and the same probably also holds for the 
CD-Bonn. If one uses the weak-coupling expression for the gap, 

A F « 2e F e~ v ^> k ^ N W , (18) 

where ep is the Fermi energy and N(0) the density of states at the Fermi level, one sees that 
the gap vanishes where the interaction goes to zero. In our phase-shift approximation, this 
happens where the phase shifts change sign, at k F « 3.6 far 1 . The Ar gonne Vis gap then 
seems to disappear somewhat too early, while the other potentials give gaps which exist up 
to what is probably unrealistically high densities. 

VI. CONCLUSION 

We have presented new calculations of the pairing gap in the 3 P 2 - 3 F 2 channel for pure 
neutron matter as a function of density. With these calculations we have aimed at establish- 
ing on a firm basis the numerical value of the gap once the bare nucleon-nucleon interaction 
is used as the pairing interaction, since in this context contradictory results have been pre- 
sented in the literature. Three different numerical methods to solve the pairing gap have 
been employed in this papery Since all three methods gave the same results, the pairing gaps 
we have obtained should be reliable from a technical point of view. 

However, our calculations have revealed that the behavior of the 3 P 2 - 3 F 2 gap at densities 
above kp ~ 2.0 fm -1 , corresponding to p ~ l-7p , where p is the nuclear matter saturation 
density, must be considered as largely unknown. Up to this point the gap is increasing (the 
values at kp = 2.0 fm -1 are about 0.6 MeV with free single-particle spectrum, and about 0.3 
MeV with BHF spectrum, independent of the potential), but how far in density this increase 
continues, depends on the individual potentials, in line with their extrapolations of the 3 P 2 
phase shift predictions. Bearing in mind that the Nijm-II potential fitted the empirical 3 P 2 
phase shift rather well up to ps 600 MeV (kp ~ 2.7 fm -1 ), we can deduce from Fig. f| 
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that the maximum gap with a free spectrum is probably below 1 MeV. How high up in 
density the gap exists must be left as an open question, although the phase shifts indicate 
that the gap should disappear at around kp = 3.6 fm _1 , corresponding to p ~ 10p - At this 
point also the purely nucleonic treatment of the dense medium is surely inappropriate. 

Before a precise calculation of the 3 P 2 - 3 F 2 pairing gap can be made, one therefore needs 
a nucleon-nucleon potential that fits the phase shifts up to m 1 GeV accurately. To us, 
the construction of potential models in which the inelasticities above = 350 MeV due 
to the opening of the NA channel are taken into account, seems to be more urgent than the 
evaluation of polarization effects on the 3 P 2 - 3 F 2 gap with the existing potential models. 
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TABLES 



k F (far 1 ) 


Bonn B 


Paris 


V U 


CD-Bonn 


v 18 


Nijm I 


Nijm II 


1.2 


0.05 


0.04 


0.05 


0.03 


0.04 


0.03 


0.03 


1.4 


0.16 


0.15 


0.19 


0.11 


0.14 


0.12 


0.12 


1.6 


0.35 


0.32 


0.45 


0.27 


0.31 


0.27 


0.27 


1.8 


0.52 


0.49 


0.75 


0.45 


0.49 


0.47 


0.45 


2.0 


0.66 


0.57 


1.02 


0.64 


0.62 


0.69 


0.68 


2.2 


0.67 


0.49 


1.14 


0.77 


0.65 


0.91 


0.90 


2.4 


0.58 


0.30 


1.13 


0.86 


0.56 


1.12 


1.15 


2.6 


0.39 


0.10 


0.95 


0.85 


0.37 


1.26 


1.39 


2.8 


0.21 




0.70 


0.78 


0.17 


1.38 


1.66 


3.0 


0.06 






0.61 


0.02 


1.37 


1.90 



TABLE I. Collection of 3 P2- 3 -^2 energy gaps (in MeV) for the various potentials considered in 
this paper. Free single-particle energies have been used. 



k F (far 1 ) 


Bonn B 


Paris 


V U 


CD-Bonn 


v 18 


Nijm I 


Nijm II 


1.2 


0.05 


0.04 


0.05 


0.04 


0.04 


0.04 


0.04 


1.4 


0.16 


0.11 


0.18 


0.10 


0.10 


0.10 


0.10 


1.6 


0.34 


0.22 


0.38 


0.18 


0.17 


0.18 


0.18 


1.8 


0.52 


0.26 


0.60 


0.25 


0.23 


0.26 


0.26 


2.0 


0.64 


0.22 


0.74 


0.29 


0.22 


0.34 


0.36 


2.2 


0.65 


0.14 


0.75 


0.29 


0.16 


0.40 


0.47 


2.4 


0.56 


0.01 


0.66 


0.27 


0.07 


0.46 


0.67 


2.6 


0.37 




0.42 


0.21 




0.47 


0.99 


2.8 


0.19 




0.23 


0.17 




0.49 


1.74 


3.0 


0.02 




0.08 


0.11 




0.43 


3.14 



TABLE II. Collection of 3 P2- 3 -^2 energy gaps (in MeV) for the various potentials considered in 
this paper. BHF single-particle energies have been used. 
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FIGURES 
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3 P2 phase shift predictions of different potentials up to E) 
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1.1 GeV, compared with 



the phase shift analysis of Arndt et al. [27]. The "old" potentials are denoted by different symbols; 
the "modern" potentials by different line styles. 
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I I I I I I 1_ 

0.5 1 1.5 2 2.5 3 3.5 




1.5 2 , 2.5 

Mfm- 1 ) 

FIG. 3. The diagonal part of the neutron- neutron potential in momentum space [Eq. @] for the 
3 i-2 (top panel) and the 3 F2 (bottom panel) partial waves obtained with the CD-Bonn, Nijmegen 
I and II and Argonne Vis potentials. 
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O 0.4 



0.2 




FIG. 4. Top panel: The angle-averaged 3 i-2- 3 -^2 gap in neutron matter depending on the Fermi 
momentum, evaluated with free single-particle spectrum and different nucleon-nucleon potentials. 
Central panel: The gap evaluated with BHF spectra. Bottom panel: The gap with the CD-Bonn 
potential in different approximation schemes. 
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FIG. 5. Effective masses derived from various interactions in the BHF approach. 
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FIG. 6. Pi gap calculated with separable potentials constructed directly from the P2 phase 



shifts. 
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